The central image of a gravitationally lensed quasar 
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A galaxy can act as a gravitational lens, producing multiple images of a background 
object. Theory predicts there should be an odd number of images^ but, paradoxi- 
cally, almost all observed lenses have 2 or 4 images. The missing image should be faint 
and appear near the galaxy's center. These "central images" have long been sought 
as probes of galactic cores too distant to resolve with ordinary observations. '^'^'^'^ '^ 
There are five candidates, but in one case the third image is not necessarily a central 
image,®'^'^'' and in the others, the central component might be a foreground source 
rather than a lensed image. ^^■^^'^'^'^^'^^ Here we report the most secure identification 
of a central image, based on radio observations of PMN J1632 0033, one of the latter 
candidates. Lens models incorporating the central image show that the mass of the 
lens galaxy's central black hole is <2xlO*^ solar masses (M©), and the galaxy's surface 
density at the location of the central image is >20,000AfQ per square parsec, in agree- 
ment with expectations based on observations of galaxies hundreds of times closer to 
the Earth. 



Central images are often called "odd" images, 
but the key property that makes them interesting 
is not their odd number, but rather their proxim- 
ity to the center of the lens galaxy. To emphasize 
the distinction, we note that there can be an odd 
number of images without any central images, if 
the gravitational field has a large quadrupole^ or 
if there are multiple lens galaxies^^'^^. 

The absence of central images is a prob- 
lem dating to the earliest days of galaxy lens 
observations, and many solutions have been 
proposed. ^'^'^'^^ Most recently, Evans &: Hunter^ 
and Keeton^ argued that the absence of central 
images at current detection limits is no longer 
surprising, given recent Hubble Space Telescope 
observations of nearby galaxies. Those observa- 
tions showed that the distribution of stars near 
the centers of massive elliptical galaxies (the pre- 
dominant type of lens galaxy) is highly concen- 
trated. Because the central image flux depends 



inversely on the square of the surface density, the 
concentrated density profiles should cause cen- 
tral images to be very faint (or even to have zero 
fiux, if the density is singular), and detection lim- 
its will need to be improved by a factor of 10-50 
before central images are commonly seen. In the 
meantime, the most favorable systems are radio- 
loud double quasars with large fiux ratios. Radio 
waves are not extinguished by dust or overpow- 
ered by light in the lens galaxy, and asymmet- 
ric doubles produce the brightest possible central 
image for a given mass distribution. 

One such system is PMN J 1632-0033.^'' It 
has two images (A and B) of a radio- loud quasar 
at redshift 3.42, with a fiux ratio of 13. The lens 
redshift is unknown; our working assumption is 
zl = 1-0, as estimated by requiring the galaxy's 
photometry and mass to be consistent with the 
fundamental plane of elliptical galaxies. There 
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is also a faint radio component (C) with a po- 
sition and flux appropriate for a central image. 
However, as in a few other systems/^ '"^^'^^ the 
possibility could not be excluded that C is an 
active galactic nucleus (AGN) in the lens galaxy 
rather than a third quasar image. The simplest 
test is to compare the continuum radio spectra 
of the components. Because lensing preserves 
the frequency of photons, lensed images have 
the same spectrum (in the absence of differen- 
tial propagation effects) whereas there is no rea- 
son why a foreground AGN would have the same 
spectrum as a background source. 

Previously we attempted this comparison at 
frequencies from 1.7 to 15 GHz, finding that C 
was fainter at low frequency than expected for 
a third image. However, the discrepancy was 
limited to a single measurement at the lowest 
frequency, where radio propagation effects are 
strongest. A central image might be affected 
more than other images by scintillation or ab- 
sorption, due to its passage through the dense 
galactic center. Thus we could make no firm con- 
clusion before obtaining data at higher frequen- 
cies, where propagation effects (scaling charac- 
teristically as z^"^) arc negligible. 

We have now extended our measurements to 
22 and 43 GHz, and obtained additional data 
at 8 and 15 GHz, using the Very Large Array 
(VLA). The high-frequency spectrum of the cen- 
tral component agrees well with those of the 
bright quasar images. For u > 1.7 GHz, the 
logarithmic slopes of flux density ratio vs. fre- 
quency (which should be zero for lensed images) 
are 0.00±0.04 for B/A and -0.02±0.07 for C/A. 
This is powerful evidence that C is a third quasar 
image. 

The evidence that C is not only a third im- 
age, but also a long-sought central image, is 
its proximity to the center of the lens galaxy 
(;S30 milli-arcseconds) , and its faintness (0.41% 
the flux density of A). This sets PMN J1632- 
0033 apart from the only other three-image sys- 
tem known, APM 08279+5255, in which the lens 
galaxy has not been detected, and the fluxes of 



all three images are of the same order of magni- 
tude. This leaves open the possibility that the 
third image in that system is not a central im- 
age, but rather a "naked-cusp" image due to a 
highly flattened mass distribution.^'^ 

For PMN J1632-0033, the central image 
properties can be used to constrain the core 
structure of the lens galaxy. In general this re- 
quires detailed modeling, for which we refer the 
reader to Ref. 15, in which some consequences 
of the central-image hypothesis were worked out 
before the status of C was clarified. Here we 
quote one result, and elaborate upon it to in- 
clude the effect of a central black hole. If the 
mass distribution is taken to be a spherical power 
law, p(r) DC r~^, embedded in an external shear 
field to account for non-sphericity, then the data 
require (3 = 1.91 ib 0.02 (2cr confidence), only 
slightly shallower than an isothermal (/3 = 2) 
mass distribution. The uncertainty is far smaller 
for this system than for typical two-image or 
four-image systems, because of the sensitive de- 
pendence of the central-image magnification on 
the exponent of the central cusp. 

If the central mass were too large, then ei- 
ther a bright fourth image would be produced, 
or the central image would be de-magnified out 
of existence.^ '^'^'^^ This fact can be used to de- 
rive an upper limit on the mass of any cen- 
tral black hole. We added a central point mass 
to the power-law model described above, find- 
ing that the limits on /3 are barely affected and 
Mbh < 2.0 X IO^Mq. 

How does this compare to the mass of the 
black hole that is expected to reside in this 
galaxy? Among nearby galaxies, Mbh is cor- 
related with the velocity dispersion of the sur- 
rounding stellar bulge. ^^'^^ The exact correlation 
is a matter of debate, but for present purposes 
we adopt the correlation of Tremaine et al?'^ 
Since the mass distribution of the lens galaxy in 
PMN J1632-0033 appears to be nearly isother- 
mal, the A-B angular separation = 1.46 arc 
seconds can be used to estimate a via the relation 
M/2 = (47raVc^)(-DLs/£'s), where £»ls and L>s 
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are the lens-source and observer-source angular 
diameter distances. Assuming a flat flu = 0.3 
cosmology with Hq = 72 km Mpc~^ and 
zl = 1.0, the results are a = 224 km s~^ and 
Mbh = (2.1 ± 0.4) X 10^ Mq. Thus, the black 
hole cannot be much larger than expected from 
the local correlation. 

A simple argument can also be used to place 
a lower bound on the core density of the galaxy. 
Since the density presumably docs not increase 
with radius, the minimum central density oc- 
curs for the case of a constant-density core, for 
which the magnification of C is //c = (1 — ^c)"^- 
Here, k is the surface density, in units of the 
cosmology- and redshift-dcpcndcnt critical sur- 
face density for lensing. The near-isothermal 
models give = 0.0082, from which it fol- 
lows KQ > 10. Under the same cosmologi- 
cal assumptions as above, the surface density is 
>20,000 Mq pc-2 within the radius of C. The ra- 
dius of C is not known accurately but is ^30 mas, 
corresponding to ^230 pc, or <0.15i?eff5 where 
i?efr is the effective radius of the galaxy. 

How does this compare to nearby galaxy 
cores observed with HST7 Among the sample 
of Faber et a/.^° are 38 early-type (E/SO) galax- 
ies with measurements of core STirfacc brightness 
profile, i?efr, and velocity dispersion. We used 
their estimates of the mass-to-light ratio (de- 
rived from the velocity dispersion) to convert the 
STirfacc brightness profile into a surface density 
profile. At R/Res = 0.15, the surface densities 
range from 10^ to 10^ MqPC"^, showing that 
our result is consistent with local estimates and 
astrophysically relevant. 

Thus we have derived limits on surface den- 
sity and black-hole mass of the lens galaxy (co- 
moving distance ^^3 gigaparsecs) that are con- 
sistent with more detailed measurements possi- 
ble only for nearby galaxies (~30 megaparsecs). 
The lensing constraints have the advantage of de- 
pending directly on mass, rather than light. This 
motivates the discovery of additional central- 
image systems, and further observations of this 
system, some of which we describe below. 



The most immediate challenges are direct 
measurement of the lens redshift and more ac- 
curate (^5 mas) measurement of the position of 
C relative to the lens center. These data would 
better pin down the physical radius within which 
the limits on surface density apply, and firm up 
the correspondence with nearby galaxies. 

Measuring correlated variability between C 
and another image would provide irrefutable ev- 
idence that C is a central image, should doubt 
still remain. Because of the different light-travel 
times, quasar flux variations should be seen first 
in A, then Ril80 days later in B, and then Ri9 
hours later in C. The ratio of the delays would 
also further constrain mass models, although the 
C-B delay measurement would be difficult due 
to the faintness of the image and the necessarily 
fine time-sampling. 

Finally, an additional image^^ is predicted to 
occur over the entire allowed range of parame- 
ters of our models with central black holes. In 
general, the fourth image has <10% the flux of 
C, although there are very rare arrangements in 
which C represents a pair of merging images. It 
is separated by <20 mas from C, indicating the 
need for very-long-baselinc interferometry. Cur- 
rently, the strongest upper limit on the flux of an 
additional component is <30% that of C, based 
on the 5(7 limit at 8 GHz.^^ Detection of the 
fourth image would be interesting because the 
separation and magnification ratio between the 
central images would provide a measurement of 
Mbh, rather than an upper limit. 
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Date 


Frequency 


A 


B 


C 


R.M.S. noi 


(UT) 


(GHz) 


(mJy) 


(mJy) 


(mJy) 


(mJy beam 


2003 Jul 23 


8.46 


202.13 


16.49 


0.84 


0.033 


2003 Aug 02 


8.46 


209.82 


16.65 


0.74 


0.040 


2003 Aug 30 


8.46 


195.19 


15.86 


0.78 


0.042 


2003 Jim 10 


14.94 


173.27 


12.50 


0.72 


0.16 


2003 Jul 01 


14.94 


180.17 


13.87 


0.57 


0.13 


2003 Jun 12 


22.46 


149.54 


10.78 


0.62 


0.076 


2003 Jun 24 


22.46 


146.95 


10.81 


0.48 


0.061 


2003 Jun 26 


22.46 


146.42 


10.78 


0.66 


0.066 


200:-! Jim 15 


4:-!.;-',4 


107.25 


7.85 


0.38 


0.077 



Table 1. Radio observations of PMN J1632-0033 with the Very Large Array. The array 

was in the A configuration. The primary flux density calibrator was 3C 286 and the phase calibrator was 
J1658+0741 (except the 8 GHz observation of 2003 Jul 23, for which it was J1651+0129). Initial calibration 
was performed with standard AIPS procedures. Further self-calibration and imaging were performed with 
Difmap. Flux densities were determined by fitting a model consisting of 3 point sources to the visibility 
data. For the; 8 GHz and 15 GHz data, the relative positions of the components were fixed at the values 
determined by previous higher resolution measurements^'*. 
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Figure 1. A radio map of the three components of gravitational lens PMN J1632-0033. The observations 
were carried out with the VLA at 43 GHz, with a total bandwidth of 100 MHz. Standard phase and amplitude 
calibration were applied, and the map was created with natural weighting. The root-mcan-squarcd (RMS) 
noise level in the residual map is <j = 0.077 mJy beam~^, as compared to the theoretical minimum RMS of 
0.04 mJy beam~^. Contour levels begin at 2.5cr and increase by powers of \/2. The synthesized beam (full 
width at half-maximum of 0.073 x 0.044 arc seconds) is illustrated in the lower left corner of the map. 
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Figure 2. The central radio component and the bright quasar images have similar radio spectra, (a) 
The flux density of each component as a function of frequency, (b) The flux density ratios relative to A. 
Results for C/A were multiplied by 10 for display purposes. All the observations at a given frequency from 
Table 1 and Ref. 15 were combined by averaging the flux density ratios and adopting the total flux density 
of the most recent measurement, to avoid problems due to source variability and inconsistencies in absolute 
flux density scale. The only discrepancy between C and A is at 1.7 GHz, where C is fainter than expected for 
a lensed image. One plausible explanation is free-free absorption due to ionized material <200 parsecs from 
the lens galaxy nucleus, the approximate position of C. The dashed line is a 2-parameter fit to the function 
/UcAe"^", using a standard approximation^^ for free-free opacity, t^, = 0.08235 T~^-^^u~'^-^E, where v is the 
frequency (in GHz) and T and E are the temperature (in Kelvins) and emission measure (the integral of 
over the path length s, in cm~^pc) of the ionized medium. The required opacity = 0.7 at the lens-frame 
frequency v « 3.4 GHz could be produced by circumnuclear material with, for example, T = 6000 and 
E = 10^ (rie = 10^ cm~^, s = 10 pc), which seems physically reasonable. For comparison, the emission 
measure is larger than in most individual Galactic H II regions (iJ~10^~^; e.g., ref. 25), comparable to some 
ionized clouds around Sgr A (i?^10®^^; e.g., ref. 26), and smaller than observed around the nuclei of some 
radio galaxies {E^IO^, e.g., refs. 27,28). The free-free absorption hypothesis and the nature of the absorbing 
medium could be tested further with future low- frequency observations. 



